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The CKM matrix element \Vus\ can be extracted from the experimental measurement of semilep- 
tonic K ^ 71 decays. The determination depends on theory input for the corresponding vector 
form factor in QCD. We present a preliminary update on our efforts to compute it in A^j = 2 + 1 
lattice QCD using domain wall fermions for several lattice spacings and with a lightest pion mass 
of about 170MeV. By using partially twisted boundary conditions we avoid systematic errors 
associated with an interpolation of the form factor in momentum- transfer, while simulated pion 
masses near the physical point reduce the systematic error due to the chiral extrapolation. 
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1. Introduction 

In the Standard Model (SM) the unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix contains 
information on the strength of flavour-changing weak decays as well as information on correlations 
between different such processes. Inconsistencies in the CKM-picture would indicate the presence 
of new physics beyond the SM. One therefore tries to determine all CKM-matrix elements as 
precisely as possible by studying flavour changing processes both experimentally (e.g. at NA62 
and LHCb at CERN) and theoretically. Here we concentrate on the determination of the matrix 
element \Vus\ from the study of semileptonic kaon (7^/3 ) decays which enables a test of first-row 
CKM-unitarity, + + = 1- is neglected because it is smaller than the errors 
in the other two terms on the l.h.s., \Vus\ and \Vud\ is known very precisely from neutron 

j8 -decay and matching the corresponding precision for is crucial in searching for deviations 
from CKM-unitarity and for possible signs of new physics. 

In this talk we present an update on our precision-study of the Ki^ form factor, /+(0) [1, 2]. 
The novelties are simulations at lighter quark masses corresponding to a lightest pion mass of about 
170MeV [3] and two additional lattice spacings. Both allow for considerably improved control over 
systematic effects. 

2. Kaon semileptonic decays 

To date, one of the most precise determinations of |Vi^^| comes from K nlv semileptonic 
decays (cf. [4]): The product |K,p|/+(0)P 

can be determined with a precision at the per mil level 
from the experimental decay rate [5] and lattice QCD provides the value for /+(0). Currently the 
level of precision set by experiment sets the precision goal for lattice simulations. 

The form factor /+ (0) is defined from the vector part of the strangeness-changing weak current 
(Yii = according to 



with /o(0) /+(0). In the SU(3) flavour limit where m\ m|, /+(0) = 1 by vector current 
conservation. It can be expanded in terms of the meson masses as /+(0) = 1 + /2 + . . . where 
fn = 0{m^^m^^m^) [6]. /2 is a known function of meson masses and the SU(3) pseudoscalar 
decay constant. In lattice simulations therefore effectively only the small higher order correction 
A/ = /+(0) — (1 +/2) is computed and extrapolated to the physical point. 

We compute the matrix element in (2.1) in terms of the ground- state contribution to suitable 
ratios of Euclidean two- and three-point functions. In a finite lattice box with periodic boundary 
conditions for the quark fields, the matrix element can in this way only be computed for meson 
momenta corresponding to the Fourier modes, i.e. with = ±0, ±1, ... for a spatial volume 
y = L^. The form factor at = is then computed by interpolating between the data for the 
form factor at these Fourier-points [1] thereby introducing a dependence of the final result on the 



(2.1) 



where 




(2.2) 
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interpolation model; different ansatze may lead to different results for the form factor, introducing 
systematic uncertanities. 



3. Partially twisted boundary condition 

In contrast to periodic boundary conditions, twisted boundary conditions allow access to 
hadron momenta other than [8]. Here we employ partially twisted boundary conditions [9] 
by applying the twist only to the valence quarks. 



(3.1) 



(6i is the twist angle in the /-direction). In our simulation a charged meson of mass m with one of 
the valence quarks twisted with angle then obeys the dispersion relation 



(3.2) 



up to exponentially supressed finite volume corrections [9]. For matrix elements like K^n, where 
the initial and final state mesons carry twists 0/ and 0/, respectively, the momentum transfer can be 
written as [10] 



= {Pi - Pff = Ei{pi) - Ef{pf) 




(3.3) 



By adjusting the twists on the initial and final meson (here tt), we can evaluate /+ directly at = 
0. The commonly used implementation is to contruct ratios of correlation functions, Rkk{pk^P7z) 
(see e.g. [10]), where either the kaon or pion is kept at rest and where the twist on the other meson 
is tuned in order to obtain results [10]: 



Rkk{pk,0) with 1 0^1 =L 
and RKni^^Pn) with \9ji\=L 



— mj. and 0;^ = , 



2mK 



(3.4) 



- m| and 0^^ = . 



The form factor /o(0) can then be evaluated directly at = as 

Rkk {pk, 0) {niK -Ej,)- Rkk (0, ^tt) {Ek nin) 



MO) 



{EK-\-mj^){mK-Ej^)-{mK + Ej^){EK-mj^) 



(3.5) 



This expression is derived by considering only the time-component of the vector current matrix 
element (2.1). Using all the other components of the weak vector current, we can construct an 
over-constrained system of equations 

RK7t0KAVt) = f+{0){EK + m„)+f-{0){EK-m„), 
RKn{0A,Vt) = f+{0){mK + E„)+f-{0)imK-E„), 

RK^tieKAVi) = f+io)dK,i+f-io)dK,i, 

RKA0A,Vi) ^ f+{0)B„,i-f-iO)9„j, {i^x,y,z), (3.6) 
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Label 


Size 


Sg 


a-' 


m;r(MeV) 


24Coarse 


24^ X 64 X 16 


Iwasaki 


1.75(4) 


330, 420, 550 , 670 


32Fine 


32^ X 64 X 16 


Iwasaki 


2.31(4) 


290, 350, 400 


32Coarse 


32^ X 64 X 32 


Iwasaki+DSDR 


1.37(1) 


170, 250 



Table 1: A summary of the three ensembles used in this analysis. Here 'Sg' denotes the Gauge action, 'm^' 
the pion mass and the lattice spacing. The 32Fine and 32Coarse data are new in this calculation. 




Figure 1: Plot of /+(0) vs. /-(O) for kaon(left) and pion(right) only twisted equations 



which we solve to obtain the form factors directly at = 0. Note that distributing the twist over all 
spatial directions maximises the number of non-trivial equations. Having eliminated the systematic 
error due to the interpolation, the remaining dominant source of systematic uncertainty is due to 
the chiral extrapolation of the lattice data to the physical point [2]. 



4. Results 



In [1] we presented results for the form factor from simulations of domain wall fermions 
(DWF) with the Iwasaki Gauge action, an inverse lattice spacing of l.TGeV and ntj^ in the range 
330-700 MeV. Here we extend the analysis by adding new ensembles 32Fine [11], 32Coarse [12] 
(cf. table 1) providing data for considerably lighter pion masses down to 170 MeV and two addi- 
tional lattice spacings. 

Closer to the physical point the twist- angles prescribed by eqs. (3.4) are very large, in particular 
for the case where the kaon is moving and the pion is at rest. This manifests itself in larger statistical 
fluctuations for the ratio Rk7i{pki 0). In order to find a better choice of kinematics we analysed the 
situation further: From eq. (3.6), the slope of /+(0) with respect to /-(O) is given by 



^/+(0) 



df-{0) 



df-{0) 



Ek + Mt, 



for the matrix element of the time-component of the vector current, \4, and 



o'/+(0) 



^/+(0) 



(4.1) 



(4.2) 
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for the spatial components V^j^z- The solutions for all equations for and (0,^;^:) on the 

32Coarse ensemble with m;^ = 250 MeV are shown in the l.h.s. and r.h.s. plots in Fig 1, respectively. 
While all solutions have a negative slope for the case where only the kaon is twisted, there are 
solutions with opposite slopes in the case where the pion is twisted. Because the solution is given 
by the intersection of the individual constraints, the kinematical situation where the pion is moving 
(twisted) and the kaon is at rest provides the best result. The statistical errors are also smaller in 
this case. Motivated by these observations, we computed all correlation functions once again for a 
third choice of kinematics with = 0, where both the kaon and the pion are twisted. This leads 
to a good constraint for /+ (0) as shown in the left plot in Fig. 2. The result is in agreement with 
the result obtained by solving all simultaneous equations for the cases where either the pion or the 
kaon are twisted as shown in the right plot of Fig. 2 (obtained by combining the plots in Fig. 1). 

We now turn our attention to the mass- and momentum-dependence of the results. In principle, 
using partially twisted boundary conditions, one is independent of the momentum interpolation. We 
wish however to include our earlier data sets, the three heavier ensembles in the 24Coarse set, for 
which we have not generated data directly at = 0. From [1] we know that in these cases the 
interpolation introduces hardly any model-dependence. In addition to the momentum-dependence 
a fit ansatz should also incorporate the SU(3) -symmetry-breaking nature of /+(0) and the strange 
quark mass dependence since the simulated strange quark mass is not exactly at the physical point. 
Our ansatz with four fit parameters Aq, Ai ,Mo,Mi [1] is 

l-q^/{Mo+Mi{mj. + ml))^ 

Fig 3 (left) summarizes the results of the analysis for all ensembles. The data points are the 
ones for the simulated, i.e. unphysical strange-quark mass. After correcting towards the physical 
strange-quark mass using the ansatz in eq. (4.3), all data points line up on the fit-curve in the r.h.s. 
plot in Fig. 3. Our preliminary result for /+(0) at the physical point is indicated by the light blue 
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Figure 3: Plot of /+(0) dependence on m\. Left : raw data points at the simulated (uncorrected) strange 
quark masses. Right : The data points are shifted to physical strange quark mass. In both plots, the curve 
uses a parameterisation for the kaon mass with strange quark held fixed at its physical value. 



square. At this early stage of the analysis we find that the statistical error at the physical point has 
been reduced with respect to our earlier results [2]. 

In our fit ansatz we use the decay constant in the chiral limit, /o, in the NLO term and have 
a form for the NNLO term which together make the ansatz consistent with the Ademollo-Gatto 
theorem and symmetry under interchange of nijj^ and mk. Since we do not know the precise value 
for /o, we repeat the global fit for /o = 100, 115 and 131 MeV. This variation in /o is also assumed 
to account partially for our NNLO term not being the full chiral perturbation theory expression. 
The result for /+(0) using /o = 115 MeV is taken as our best value and the values at 100 and 
131 MeV provide an estimate of the systematic error from the chiral extrapolation. 



f+(0) 



RBC/UKQCD 12 PRELIMINARY 
JLQCD 11 
RBC/UKQCD 10 
RBC/UKQCD 07 

FLAG-1 estimate for Nf =2 + 1 Eur. Phys. J. C (2011) 



ETM lOD 
ETM 09A 
QCDSF 07 
RBC 06 
JLQCD 05 
JLQCD 05 

FLAG-1 estimate for Nf =2 Eur. Pliys.J. C (2011) 



0.94 0.96 0.98 1.00 



Figure 4: Comparison of recent Lattice results 
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In Fig. 4 we compare our new preliminary result for /+ (0) with other recent Lattice determina- 
tions. We emphasize that our result is preliminary, but we expect to reach a new level of precision 
when the analysis is complete. 

5. Conclusion 

We have presented an update on our simulations aiming at the precise prediction of the kaon 
semileptonic form factor at vanishing momentum transfer. The novelties are simulations for pion 
masses down to 170MeV and for three different lattice spacings. We identified preferred choices 
for the kinematics when using partially twisted boundary conditions in order to simulate directly 
at = 0. In particular, we now understand why twisting either only the pion or both the kaon and 
the pion can constrain the form factor better than twisting only the kaon. At this early stage of the 
analysis the inclusion of the new ensembles in our global fit leads to a reduction of the statistical 
error at the physcial point compared to our earlier result [2]. This indicates significant progress in 
the study of K ^ tt semileptonic form factors which would directly impact our knowledge of the 
CKM Matrix and provide for improved constraints for new physics beyond the Standard Model. 
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